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ABSTRACT 



Context. The carbon-to-oxygen (C/O) ratio of asymptotic giant branch (AGB) stars constitutes an important index of evolutionary 
and environment/metallicity factor. 
' Aims. We develop a method for mass C/O classification of AGBs in photometric surveys without using periods. 

Methods. For this purpose we rely on the slopes in the tracks of individual stars in the colour-magnitude diagram. 
, Results. We demonstrate that our method enables the separation of C-rich and O-rich AGB stars with little confusion. For the 

Magellanic Clouds we demonstrate that this method works for several photometric surveys and filter combinations. As we rely on 
no period identification, our results are relatively insensitive to the phase coverage, aliasing, and time-sampling problems that plague 
period analyses. For a subsample of our stars, we verify our C/O classification against published C/O catalogues. With our method 
we are able to produce C/O maps of the entire Magellanic Clouds. 

Conclusions. Our purely photometric method for classification of C- and O-rich AGBs constitutes a method of choice for large, 
near-infrared photometric surveys. Because our method depends on the slope of colour-magnitude variation but not on magnitude 
zero point, it remains applicable to objects with unknown distances. 

Key words. Methods: data analysis, Techniques: photometric, Surveys, Stars: oscillations (including pulsations), Stars: AGB and 
^— ^ post- AGB, Magellanic Clouds 
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C> 1. Introduction It is believed that the rate of conversion of O-stars into 

' C-stars depends on the efficiency of the third dredge-up and 

c ^ Low and intermediate mass stars evolve through three late me extent an d time variation in the massloss (e . g . Iben 1981; 

^ , stages. After passing the red giant branch (RGB) they reach Marigo et al . 1999) . Mass loss is expected to be stronger in 

^ ■ maximum luminosity at its tip (TRGB) followed by a luminos- metal . rich starSi ivi t0 shorter AGB life and eventually yield- 

L • ; ity drop after the helium flash. Next they grow again and move ■ a c . star Thefje thermal lfjes resuk fa ^ luminosit varia . 

along the asymptotic giant branch (AGB). According to their tionfj wkh the peak . t0 .peak amplitude up to a few magnitudes at 

spectra, the AGB stars split into oxygen-rich stars (O-rich, M- v j sua j wave lengths 

stars, or K-stars) and carbon-rich stars (C-rich, C-stars, or N- ' 

C3 ■ stars). The M-stars have more oxygen than carbon in their at- According to llben & Renzinil (ll983l). the lower the metallic - 

mospheres (C/O < 1). When the abundance of oxygen equals ity the less carbon needed to be dredged-up in order to convert 

that of carbon the AGB type is S. Classification of the lat- an O-star into a C-star, hence the correlation between metallicity 

ter is difficult and also involves intermediate MS and SC types and the C/M ratio. For lower metallicity, the AGB evolutionary 

dCioni & H abing 20031). tracks move to higher temperatures; for very low metallicity a 

Stars entering the AGB phase are rich in oxygen, how- post-horizontal-branch star may become a white dwarf without 

ever, subsequent dredge-up caused by thermal pulsation may first becoming an AGB star. Since O and C stars reveal very dif- 

enrich their atmospheres in carbon (llben & Renzinilll983l) . On ferent spectra, it is relatively easy to spectroscopically identify 

the AGB, several thermal pulses (TP- AGB) may occur, resulting C-rich stars even at large distances. The differences in molecu- 

in dredging up the matter enriched in carbon nuclei by nuclear lar blanketing and dust creation result in an observed sharp di- 

fusion. In this way the O-rich stars are converted after at least chotomy in infrared colours of O- and C-rich stars. Thus the O- 

several pulses into the C-rich ones (iMarigo et al.ll2008l and ref- and C-stars are often identified using narrow filters. The J - K 

erences there). colour of the C-rich stars is > 1 .4 mag and systematically red- 

der than that of O-rich s tars (lFrogelet al.1119901 ICosta & Frogell 

Send offprint requests to: A. Schwarzenberg-Czerny fl996tlCioni etal , 2000a). These effects influence the luminosity 
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function of O- and C-rich AGB stars. The Magellanic Clouds are 
good test-beds of theories of the late stages of stellar evolution. 
They are nearby and yet far enough to ignore the LMC thickness, 
so that all stars are approximately at the same distance. 

Late evolutionary stages are often associated with long 
period- (LPV), semi-, and irregular-variability. A wealth of 
data on this variability was collected as a by-product of 
the microlensing surveys. All yielded catalogues of numer- 
ous variable stars. In these data the period-luminosity rela- 
tions of the LPV's on the AGB are well documented for both 



the L arge a nd Small Magella nic Clouds (Woodetal. 199' 
modi [20001: ICioni et all [20011 |20(jl iNoda et alj 120021 12 004 
Lebzelter & Hinkld 120021: llta et all l2004allbfc iKiss & Beddin 



2003L l2004t|Soszynski et al.i2004al|bl l2005l:[Groeneweger]200 



Fraser et al. 2005; Raimondo et al. 2005) 



The period-luminosity diagram for red variables at advanced 
evolutionar y stages reveals six sequences for different classes 
of objects dlta et alj 12004b). Miras and some low-amplitude 
semiregulars pulsating in the fundamental mode form the se- 
quence C. Other semiregular stars pulsating in the second and 
third overtone modes occupy sequences A, B, and C, and the 
RGB eclipsing binaries form their own sequence E. The origin of 
the sequence D corr esponding to long- secondary periods (LSP) 
remains ambiguous. IWood etal.l d2004) considered a number of 
possible explanations of LSP (radial and non-radial pulsations, 
rotations, orbiting companions, chromospheric activity, orbiting 
dust clou ds), but none fit th e observations satisfactorily. 

ISoszvnski etal] d2004al) combined the OGLE-II and OGLE- 
III data and found the multiperiodicity of the red giants to be 
variable with a small amplitude. They exhibited two modes 
closely spaced in their power spectrum. This is likely to indi- 
cate non-radial oscillations. They also show that members of the 
short-period P-L sequences below the TRGB c onstitute a mix- 
ture of the RGB and AGB varia bles. Recently. ISoszvnski et ail 
(2004a, 2005l): ISoszv6skil (120071) . and lDerekas et al.l(l2006l) have 
demonstrated that the sequence E overlapped with the sequence 
D, which may be evidence for a binary origin of the sequence D. 

There are differenc es between the O-rich and C-rich LPVs, 
too. ICioni et alJ (120031) note th at C-stars have a larger amplitude 
than O-stars. llta et al.1 d2004al) confirm that O- and C -rich Miras 
follow different period vs. (J - K) colour relations dFeast et al.l 
1989). The I-band amplitude of C-rich Miras tends to grow with 
the redder mean ( J - K) col our, while the amplitude of O-rich 
Miras is colour independent (Ma tsunaga et al.l l2005). 

The evidence has provided new and significant con- 
straints for theoretical pulsation models. Models and obser- 
vatio ns are compared in three ways: using stellar isochrones 
(e.g. iBressan et all 119961: iBruzual & Chariot! 120031: iMouhcind 
2002; Marig o etal.1 120031: IGroenewegen & de Jongl Il993l) 
stella r tracks (IGro enewege n & de Jong _ 19931 Marigo et al.l 
1999t iMouhcine & Lanconl 120021: lLancon & Mouhcinel 120021: 



Mouhcine & Lancon 2003), and the fuel consumption theorem 
bv lRenzini & Buzzonil d 1986h and lMarastonl dl998ll2005b . 

The synthetic characteristics of AGB stars are derived either 
from complete evolution tracks, semi-empirical fits to the core 
mass-luminosity relation, or from the core mass-interpulse pe- 
riod relation and the mass-loss rate as a function of stellar pa- 
rameters. The examples of suc h c alculations are prov ided by 
IGroenewegen & de Jong] d 1993b and llzzard etal] (120041) . 

Recent attempts have been made to include the TP-AGB 
phase in evolutionary population synthesis models. Both im- 
provements in the low-temperature opacities and peculiarities 
of the actual surface chemical composition have such profound 
consequences that the whole AGB evolutionary scenario be- 



came significantly affected (Marigo etal. 2008). Use of the 
molecular opacities reflecting the actual chemical composition 
leads to a significa nt decrease in T e ff as soon as C/O > 1 
(Marigo et al. 2008 and references therein). Thi s decrease is 
confirmed by observation s of galactic AGB stars (Bergeat et al. 
l200H lMarigo et al.ll2003l) . and it naturally explains the presence 
of a red tail in the (J - K) colour-magnitude diagram of C stars. 

The mass loss is driven by pulsation in a complicated way. 
The pulsation pushes a fraction of the atmosphere above the 
photosphere, creating a cool and dense environment where dust 
grains form and grow efficiently. The radiation pressure on the 
dust combined with the momentum c arried in the shock waves 
drives the mass lo ss dWoodl Il979t iBowen & Willsonl 119911: 
iHoefner et al]|1996l e.g .) 

Lebzelter & Wood (2005) compared the observations of 



LPV in 47 Tuc with the models. On one hand, they find that 
the models without mass loss fail to reproduce the observed pe- 
riods of the small amplitude pulsators. On the other, the K-log P 
sequence of the large amplitude variables, such as the Miras, is 
inconsistent with the mass-oss models and consistent with the no 
mass-loss models. Only models involving the non-linear funda- 
mental mode yield per iods consistent with the Miras in 47 Tuc 
(lOlivier & Woodll2005l) . 

The past decade has b rought results of exte nsive pho- 
tometric surveys - OGLE (Pacz vnski etall 1 19941) . OGLE II 
([Udalskiet al l 1 19971) . EROS dAubourgetal.il 19951) . MACHO 
(lAlcocketal.1 11997). and MOA (Bond etalMOOll) - covering 
several optical and infrared bands, sometimes simultaneously. 
They have provided a wealth of observations of red variable 
stars, thus enabling the study of their population properties. In 
particular, use of infrared colours and/or pulsation periods have 
enabled the classification of C- and O-rich stars. We discussed 
above the relevant observational and theoretical results. 

In the present paper, we attempt to analyze the photomet- 
ric properties of red variables in the visual region without re- 
course to their periods. Our analysis should be complementary 
to any traditional methods while suffering less, if at all, from 
aliasing and seasonal interference. We employ correla tion slopes 
of col our and magnitude variability introduced by IWood et ail 
(12004 . In this paper the slope or of the correlation of R magni- 
tude and V - R colour variations revealed no relation with any 
other property of red variables. We employ a different filter com- 
bination in our attempt of photometric classification of C and O- 
rich red variable stars. Our results should not depend much on 
the number of pulsation cycles covered by a given survey. Our 
practical requirement that data span at least one pulsation cycle, 
is met in most large surveys. 

In Sect. [2] we describe data employed in the present study. 
Our calculation methods are described in Sect. [3] In Sect. |4]we 
propose new methods of selection O- and C-rich stars. Due to 
large differences in filters and time span, we present our results 
for each survey separately. We give an estimate of the C/O ratio 
for variable stars in the Magellanic Clouds. In Sect. [5]we employ 
a simple model to demonstrate how observed effects may arise. 
Possible peculiar variables are discussed in Sect. [6] We conclude 
in Sect. ~] 

2. The data 

Our paper is based primarily on the EROS-2 photometric survey 
of the Magellanic Clouds cross-referenced with the 2MASS in- 
frared magnitudes (Sect. |2~TT >. We explicitly indicate when these 
are supplemented by the OGLE and MACHO photometry and 
four catalogues of C-and O-rich stars (Sect. l2~2l - I2.4t . 
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2.1. EROS-2 survey 

The Experience de Recherche d'Objets Sombres (EROS-2) 
project employed extensive photometry obtained with the 1- 
meter MARLY telescope at La Silla Observatory, Chile, to 
search for the baryonic dark matter of the Galactic halo 
by means of the gravita t ional microlensing (lAfonso et all 
20031; lLasserre et all 120001; iPalanque-Delabrouille et al.l 11998; 
Tisserand et al.1 120071) . The observations were performed be- 
tween July 1996 and February 2003 through a dichroic plate- 
splitting light into two wide field cameras covering 0.7 x 1 .4° in 
right ascension and declination each, yielding two broad pass- 
bands. The so-called blue channel (420 - 720 nm, hereafter Be 
band) overlapped the standard V and R standard bands, while 
the red one (620 - 920 nm, hereafter Re band) roughly matched 
the standard I filter. Each camera constituted a mosaic of eight 
2k X 2k CCDs with a pixel size on the sky of (0.6") 2 - 

Ten fields covered the SMC and 88 fields covered the LMC. 
The photometry from individual images was combined into light 
curves using the P eida package developed specifically for the 
EROS experiment (lAnsaril 19961). The es t imated accuracy of this 
photometry is discussed by iDerue et al.l (120021) . For uniformity, 
the SMC data were analysed again with the more recent version 
of Peida. 

The RGB and AGB variable stars were selected from the 
SMC and LMC lists of EROS variables using their location in 
the colour-magnitude diagram. In this way we selected 28914 
stars in LMC and 5930 in SMC. They constitute our sample used 
in further analysis. 

2.2. OGLE data 

The Optical Gravitational Lensing Experiment (OGLE-II) data 
were collected with the 1.3 m Warsaw Telescope at the 
Las Campanas Observatory, Chile, operated by the Carnegie 
Institution of Washington. The /-band data span about 3000 
days: from January 1997 to April 2005. The V-band mea sure- 
ments were obtained from 1997 to 2001 (IZebrun et al.ll200lT) . so 
they span a shorter time baseline. Up to 70 V-band points per 
star are available. In the /-band, 5 00 to 900 measuremen ts were 
available, depending on the field (Soszvnski et al. 2005). In the 
pr esent work we u s e dat a on 3586 LPV from LMC published 
by ISoszvnski et al. (2005), downloaded from the OGLE home- 
pagtQ 

We cross-identified EROS and OGLE objects within a 1 .3" 
radius. The OGLE-II fields cover only the bar of LMC, i.e. a 
much smaller field than for EROS, hence there are relatively few 
stars in common. 



2.3. MACHO data 

The MAssive Com pact Halo Objects (MACHO) project 
dAlcocket all 1 19971) comprises eight years of observations of 
the Large and Small Magellanic Clouds and of the Galactic 
Bulge. For the present purposes we s elected LMC stars f rom the 
MACHO catalogue of variable stars dAlcock et al.ll2003l) . Using 
the MACHC0 interface, we downloaded 2868 individual light 
curves of the stars classified as the red giant variables (classified 
as Wood A, B, C, and D sequences). 



2.4. 2MASS survey data. 

The Two Micron All Sky Survey (2M AS Sf| employed two 1.3- 
m robotic telescopes, one at Mt. Hopkins, Arizona, and one at 
CTIO, Chile. Each telescope was equipped with a three-channel 
camera, capable of observing the sky simultaneously at J (1.25 
microns), H (1.65 microns), and K s (2.17 microns). The south- 
ern facility began collecting Survey data in March 1998 and con- 
ducted its final scan in February 2001. 

We extracted 2MASS data for each of the EROS fields sep- 
arately. Then cross-identification was done separately for each 
field. During the cross-identification of EROS and 2MASS, we 
employed a 3" search radius. We selected only stars with pho- 
tometry available in all three bands J,H, and K s . AGB stars are 
among the brightest stars of the LMC. If an area of the order of 
16 square degrees contains 25000 AGB stars, their average sep- 
aration is 1.5 minutes of arc, so the probability of blending two 
of them within the search radius is (1 .5 * 60/3)~ 2 = 0.0006, thus 
corresponding to fewer than 15 stars in the whole sample. This is 
negligible compared to the blur of our photometric plots. At this 
point it is worth discussing the contamination of our AGB sam- 
ple of IR colours with mistaken cross-identification with RGB 
stars. For stars brighter than TRGB (K < 12 for LMC and 
K < 12.7 for SMC), this problem does not exist. Suppose, for 
an EROS AGB star, that IR colours for an RGB one were fit- 
ted, then the star would land leftwards of the vertical lines on 
ay versus K diagrams. One possible problem is mixing faint 
AGBs with RGBs, causing pollution of RGBs with AGBs but not 
otherwise [i.e. faint AG B remain clean, if possibly incomplete]. 
iKiss & Bedding <2004l) consider the same problem in more de- 
tail for stars from OGLE and 2MASS and conclude that con- 
tamination is insignificant, possibly no more than 0.3%, for their 
search radius of 1", corresponding to less than 3% contamination 
for our radius of 3". Since there were selected MACHO objects 
than for EROS, we cross-identified MACHO and 2MASS with- 
out subdivision into fields and with the same 3" radius. Again 
we only selected stars with photometry in all three bands J,H, 
and K s . 



2.5. The catalogues of C-rich stars 

We selected 1707 C-ri ch stars from the SMC identified 
by iRebeirot et al.l d 19931) in the low-resolution spectroscopic 
survey employing the ESO 3.6 m telescope and 1185 
stars identified in the Siding S pring Observatory survey by 
iMorgan & Hatzidimitriould 19951) . 

A catalogue of 7760 C-rich stars in the LMC was presented 
by iKontizas et alj {2001). These stars were identified during a 
systematic survey of the objective-prism plates taken with the 
UK 1.2 m Sch midt Telescope. To th ese we added the list of C- 
rich stars from iGroenew egen (2004|)- 

Sosz vnski et al.l d2005l) demonstrate a new method of dis- 
tinguishing between O-rich and C-rich Miras, SRVs and stars 
with long secondary periods, relying on their V and /-band pho- 
tometry and periods. The list of stars selected from the LMC in 
this way was downloaded from the OGLE homepage. All C-rich 
stars from these catalogues were cross-identified with the EROS 
and MACHO stars within a 3" radius. 



1 http://sirius.astrouw.edu.p1/~ ogle/ 

2 http://wwwmacho.mcmaster.ca/ 



3 http ://w w w. ipac . caltech. edu/2mass/ 
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3. Data pre-processing 



The analysis outlined in this paper was performed separately for 
EROS LMC and SMC, OGLE LMC and MACHO LMC data. 
We used data from all surveys to study a path in the colour- 
magnitude diagram traversed by each variable star. We investi- 
gated correlation of the average slope of the path with the chem- 
ical composition of the objects. The results of this analysis are 
summarised in diagrams described in the Sect. |4] 

The EROS survey produced simultaneous light curves in two 
spectral bands. By applying time filters to each band light curve 
we obtained a smooth, low-pass filtered light curve and its com- 
plementary high-pass light curve. The sum of the two reproduces 
the original light curve. In this way for each spectral band we 
obtained 3 light curves: raw, low-, and high-pass, six in total. In 
further analysis we study correlations of stellar properties with 
the properties derived from these light curves. 

3.1. Photometric calibrations. 

Both EROS and MACHO obtained simultaneous expositions in 
two bands termed "red" (Re and r, respectively) and "blue" (Ve 
and v, respectively) for nearly all observations. We simply ig- 
nored observat ions made in just one filter. 

Following Tiss erand et al.l(l2007l) . we converted the raw pho- 
tometry into the standard Kron-Cousins system using the follow- 
ing relations: 



Veros 
Ieros 



1.666V £ -0.666/? £ 
Re- 



(1) 
(2) 



The MACHO data from the web database are availa ble in 
the raw instrumental system. Following lAlcock et al.l (1 19971) . we 
converted the raw photometry into the standard Kron-Cousins 
system using the following relations: 



Vmacho = v + 23.699-0.1804(v-r) 
Rmacho = r + 23.412 -0.1825(v-r) 



(3) 
(4) 



where v and r are instrumental magnitudes and Vmacho an d 
Rmacho are m the Kron-Cousins standard system. These cali- 
bration formulae are estimated to have an overall absolute ac- 
curacy of +0.1 mag in Vmacho or Rmacho an d ±0.04 mag in 
(V - R)macho dAlard et al.l l200l1l. 

The OGLE survey employs one camera, so observations in 
different filters are not simultaneous. Observations in V filter 
were obtained much less frequently than in /. For each V ob- 
servation, we searched the nearest / observations. If two / points 
before and after a V observation separated by no more than 2 
days, we performed a linear interpolation at the time of V obser- 
vation. For just one / point nearby, no more than one day from 
the V observation, we assume they were measured nearly simul- 
taneously. In this way for most V observations we found its cor- 
responding / magnitudes. Other V observations were ignored. 
Because we only analyse stars with long periods, this procedure 
proved reliable for our purposes. 

3.2. Time filtering 

Through the whole paper except for Sect. |6] we used the raw 
light curves. Light curves for Sect. |6]were filtered in the time 
domain. The low-pass curve was obtained from the raw one by 
means of the "moving median". The purpose of time filtering is 
to separate short period pulsations in some stars from the LSR 
Since our procedure does not depend on actual detection of these 




Fig. 2. Sample plot of raw colour and magnitude variations for 
the same 2 stars as in Fig. Q] We also plot the regression lines. 
Their inclination parameter is ay, discussed in the text. 



periods, we applied it to all EROS red variables. We selected 
a window of width w centred on a given point and replaced its 
value with the window median value. In practice we filtered each 
light curve twice consecutively using windows of width w\ and 
W2 dependent on average K s luminosity: 



m = 10 (*-Z»/(-3.9) 
W2 = 10 C*-^D/(-3.9) 



(5) 

(6) 



where the D magnitude is equal to 19.4 for the LMC and to 2 1 . 1 
for the SMC. These windows were selected to have intermediate 
length between the short period and LSP (see Sect.[6]for discus- 
sion of these periods). As both periods depend on mean stellar 
brightness K, so do our window widths. The slope of -3.9 cor- 
responds to a border line between the short period sequences 
and the L SP sequence in the perio d-luminosity (P - L) diagram 
(llta et al.ll2004bt ISoszvnski et al.ll2004ah . As this area is virtu- 
ally devoid of stars, its exact value is of no consequence for the 
current considerations. 

We found that a light curve smoothed with two moving 
medians reveal LSP more clearly than when smoothed only 
once. Any short-period variations are filtered out. The high pass 
curve was obtained by subtracting the low pass one from the 
raw data. Only short time variations are left in the high-pass 
light curve. Sample light curve s are plotted in Fig. Q] A simi- 
lar procedure was employed by IWood et aD (120041) except that 
we did not use the short period as the parameter determining fil- 
ter wi ndow to avoid confusio n for multi-periodic/irregular red 
giants (ISoszvnski et al.ll2004ah . For EROS and MACHO colours 
(V - I)eros , (V - R)macho were calculated for each time point. 
Subsequently we filtered colours similarly to how magnitudes 
are filtered. In this way we got the colour variation for raw data, 
LSP, and short time scale variations separately. For OGLE the 
number of V observations was too small to filter. 



3.3. Analysis of colour variations 

In further analyses we considered variations in colour against 
magnitude, for all possible combinations of light curves (raw, 
long-, and short-time scale colours against raw, long- and short- 
time scale V and / magnitudes). Sample plots are presented in 
Fig- El Next we fitted by the least squares the regression lines, 



(V-I) = ail + bj 
(V-I) = a v V + b v 



(7) 
(8) 



where ay and a/ denote the slope of the corresponding colour 
- magnitude relation. This kind of linear equation was fitted for 
EROS, OGLE and MACHO data. 

To verify the quality of our slope parameter, we calculated 
the correlation coefficient p and the statistic t in the following 
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Fig. 1. Sample EROS light curves for two stars with spectral types assigned bv lCioni etail (1200 ll) : DCMC J052618. 12-694100.0 
of type C (top row) and DCMC J05 2446.9 1-694949.9 of type M (bottom row). The left, middle, and right panels present the raw, 
low -pass, and high-pass light curves, respectively. For details of filtering see Sect. 13.21 



way: 

= c y { y,y-/) 

CTvCTv-l 

t = ■ 9 yJn- 7. where (10) 

i " _ 

Cov{V,V-I) = - Y(V- V)[V-/-(V-/)]. (11) 

n t-r 1 

i 

Results are shown in Fig [3] For the null hypothesis Ho assuming 
no correlation, i.e. p = 0, a t statistic obeys the Student distribu- 
tion with n — 2 deg rees of freedom, where « is number of points 
per star dFiszll 19631) . For a probability of 0.995, t is equal to 2.6. 
For most stars the correlation coefficient is greater than 0.5, and 
t is as significant as 15 or more. 

3.4. Analysis of amplitude variations 

The amplitude of the luminosity variation was calculated from 
the distance of the maximum and minimum in the correspond- 
ing, possibly filtered light curve. We derived the amplitude of 
the colour variation as the distance between points on these lines 
corresponding to extreme magnitudes. In this way we dimin- 
ished the influence of individual colour errors on the colour am- 
plitude. These amplitudes were calculated for each filter, and raw 
, long-, and short time-scale data. 

4. Photometric chemical classification of red 
variable stars 

4.1. The slope-amplitude diagram. 

Wood et al.l (|2004) have introduced the variability slope parame- 
ter flR. They note it has no obvious relation to any other property 
of red variable stars. Noticing that red giants have greater am- 
plitude in the blue filter results in better accuracy for the slope 
estimates. We adopted the ay slope parameter as a better substi- 
tute. In the left panels of fig. [4] we plot the slope ay of the colour 



- blue magnitude relation against amplitude, for all EROS red 
variable stars in the LMC and SMC. For completness we also 
plot ay against the near infrared magnitude and colour in the 
middle and right panels. In the rest of the present section, we in- 
vestigate our slope parameter ay as a tool for the chemical and/or 
evolutionary status classification of red variable stars. We defer a 
discussion of the classification corresponding to the middle and 
right panels until Sects. I4~3l and l4.4l Anticipating our result, we 
adopted in Fig. [4] different colour coding for O-rich RGB stars, 
O-rich, and C-rich AGB stars as green(light grey), blue(black), 
and red(dark grey) dots, respectively. The colour-coded selec- 
tion, for all diagrams, comes from the third method presented in 
the right panels. 

The left panels in Fig|4]display plots of the slope ay against 
the colour amplitude ampy; for EROS LMC and SMC data, re- 
spectively. Note the conspicious two-modal distribution of aV 
suggesting separate clustering of C- and O-rich AGB stars. In 
these plots there are two distinct branches corresponding to 
steep and shallow slopes, i.e. for large and small ay. For both 
branches, the slope initially grows and then, for amplitudes over 
0.3, saturates at the constant values of 0.65 and 0.3. These two 
groups can be separated by the curve with the equation 

fly = 0.52tanh(am/W0.15). (12) 

The same border holds for both EROS LMC and SMC data. In 
the following sections we demonstrate by cross-checking in the 
catalogues that the stars above and below the curve are respec- 
tively O- and C-rich. Thus our Eq. (fT2l constitutes the new clas- 
sification criterion for C- and O-rich stars. Its importance and 
novelty stem from using the visual photometry alone with no 
recourse to periods. Our criterion does work for an incomplete 
phase coverage and even for data spanning as small an interval 
as a typical Mira period. Moreover, at this point there are no ob- 
stacles to test the usefulness of our new method for the semi- and 
irregular red variables. Problems in using aj for low-amplitude 
variables are illustrated in Fig. [5] The two clusters merge for low 
amplitudes. In this way we confirm that the slope ay is indeed 
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Fig. 4. Plots depicting three methods of distinguishing O- and C-rich variable stars (from left to right): slope-colour amplitude 
av(ampvi), slope-K magnitude ay(K) and slope-/ - K colour ay(J - K), for EROS LMC (top) and SMC (bottom) data. Vertical 
lines indicate adopted TRGB dividing RGB and AGB stars. See text for a detailed description of the methods and parameters of the 
boundaries. Black, dark grey, and light grey dots (blue, red, and green in electronic edition) respectively mark O- and C-rich AGB 
stars and O-rich RGB ones, classified using K magnitudes for AGB/RGB and leftmost panels for O/C content. 
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much more suitable for classi fication purposes than our a/ and 
a R used by Wo od et al.1 d2004l) . 

In Sects. [4721 and !4.61 we verify our classification against ex- 
isting catalogues, while in Sects. 1431 - 14. 5l we discuss alternative 
photometric systems. 

4.2. Photometric properties of C-rich and O-rich stars. 

Of the 50 stars assigned bv lCioni et al.l (1200 ll) spectral types C, 
M, and S, 30 lie in our fields and have entries in 2MASS. In Fig. 
[6] we plot them in the same way as in Fig [4] Spectroscopic and 
photometric grouping is consistent in general and that in partic- 
ular the S-stars lie on the border between C and M stars, as ex- 
pected. Clear separation of C and M stars in the left panels con- 
stitutes the most reliable test of our method. More extensive tests 
follow from cross-identifying of the stars from EROS, OGLE, 
an d MACHO with the s pectrosco pic catalogues o f C-rich stars 
bv lKontizas et alj d2001l) andlGroenewegenl d2004h for the LMC 
and by [Rebeirot e t al.l d 19931) and iMorgan & Hatzidim itrioul 



d!995l) for the SMC (Sect. |23J. For this purpose the optical ( V, 
/) magnitudes were obtained from EROS, OGLE, and MACHO 
photometry, while infrared (/, K) magnitudes were obtained by 
cross-reference with 2MASS (Sects. O O O and H3J. In 
Fig. [7] we plot all identified C-rich stars in tha same fashion as 
in Fig. |U The border lines in Fig. [7] are plotted for correspon- 
dence with Fig. |4] All previous succesful photometric criteria 
wer e based on the infrared colours (c.f. middle and right panels, 
e.g. lSoszvflskiet al.1 d2009l) ). 

Some confusion in the distribution may be caused by pres- 
ence of the intermediate S type stars with the number of O-atoms 
equal to that of C-atoms. This confusion arises from the diffi- 
culty of classifying the S stars using the low-resolution spectra 
obtained for the catalogues employed in this section. 

4.3. The slope-luminosity diagram. 

Infrared colour and luminosity may be employed to addition- 
ally split red stars according to their evolutionary status . In the 
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middle panels in Fig. [4] we plot slope ay against K luminos- 
ity for EROS LMC and SMC, respectively. The separation of 
the C- and O-star sequences is as clearly visible as before, yet 
the O-sequence terminates abruptly near the K magnitude of 
the TRGB of 12 an d 12.7 for the LMC and SMC, respectively 
dCioni et al. l l2000cT) . The difference in magnitudes corresponds 
to the difference in distance moduli of the LMC and SMC. The 
horizontal line at ay=0.52 corresponds to the criterion of Sect. 
14.11 The manifestation of the TRGB here clarifies somewhat the 
behaviour observed in Sect. 14.11 For the O-rich RGB stars be- 
low the TRGB the slope ay grows with K magnitude. As soon 
as they reach the AGB, the slope saturates. The skew separation 
lines plotted in the figures correspond to the equations 

a v < -0.4^ + 5.31 for ^ > 12.0 (13) 
a v < 0.52 for K < 12.0 (14) 

for the LMC and 

a v < -OAK + 5.51 for K > 12.7 (15) 
a v < 0.52 for K < 12.7 (16) 

for the SMC. 

In the plots discussed in the present section, stars group in 
three clusters corresponding to the RGB and O- or C- AGB, but 



the K luminosity of the AGB and RGB stars overlaps to some 
extent. A fraction of the AGB stars exists with luminosity be- 
low that of the TRGB. Since the RGB and O-rich AGB stars 
have similar slopes, the RGB region is tainted with some O-rich 
AGB stars. In the above, we only employ K magnitudes for mor- 
phological purposes to separate stars on the colour diagrams. For 
stars with dust shells, the K magnitude fails as a bolometric lumi- 
nosity indicator; however, in the present paper we concentrate on 
optical photometry, and for the detailed discussion of IR proper- 
ties the reader is referred to the original papers. 

4.4. The slope- J - K colour diagram. 

In the infrared, the C-rich star s appear systematically redder (J— 
K > 1.4) than O-rich stars (iFroeel et al.lll990t ICosta & Frogell 
1996;; ICioni et alJl2000at ISoszvnski et aW OOS): however, in K 
vs. J - K diagrams (Fig. [8} there is some overlap between the C 
and O regions. The efficient way to separate the C and O stars is 
to combine the slope a y with the J — K colour. In the right panels 
of Fig.|4]we plot slope ay against J — K. Now separate clustering 
of C- and O-rich stars becomes more obvious. Based on this we 
propose the following criterion for C-rich stars: 

a v < 1/3/T+ 1/15 and J-K>\.\5. (17) 
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Fig. 3. Plot of the colour-magnitude correlation coefficient 
against the slope of magnitude-colour relation, for EROS LMC 
and SMC data. The correlation coefficient serves as a measure of 
quality of the slope a v . The horizontal lines correspond to prob- 
ability 0.995 of the null hypothesis Hq, i.e. points lying above 
correspond to statistically significant correlation. 



This criterion corresponds to the colour coding of stars in all 
panels in Fig. [4] Oxygen stars were separated from RGB and 
AGB using TRGB. 

The Z-like shape of the distribution of stars in the ay - (J-K) 
panels in Fig.|4]represent the evolutionary sequence from O-rich 
RGB (top left) through O-rich AGB (right) down to C-rich AGB 
objects. For O-rich RGB stars, amplitudes tend to be low (left 
panels in Fig. [U and slope parameter ay is poorly constrained. 
For higher luminosity O stars beyond TRGB, the slopes are 
better determined and stars concentrate around ay = 0.65 and 
J - K — 1.2. When TP- AGB begins and C-rich matter shows 
on the surface, the star migrates towards the inclined border line 
where there are stars with intermediate MS, S, and SC types. 
Eventually after crossing the border, it lands in the C-rich AGB 
star cluster. There is no appreciable difference between diagrams 
for the LMC and SMC. Using ay, K, and J - K we cannot only 
classify stars but also we could roughly determine its evolution- 
ary stage. No such direct evolutionary interpretation was avail- 
able for our slope parameter diagrams; however, this was made 
by resorting to cross identification of stars in the infrared and 
optical catalogues, and looking for which slope diagram con- 
centrations correspond to the concentrations in the IR diagram. 
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Fig. 8. K - J-K diagram for LMC and SMC red variables. 



4.5. The Wesenheit index-luminosity relation 

The Wesenheit index is a reddening-free parameter defined by 
a linear combination of stellar magnitudes. For EROS bands 
it coul d be calculated from the following equation (e.g. iTanvirl 
(fl997l) and references there): 



Wieros - Ieros - l-55(VgsQ5 - Ieros) 



(18) 



where Ieros an d Veros are light curve average magnitudes 
(Sect. 13. 11 1. For EROS, observations in both filters are simulta- 
neous so time sampling yields no big uncertainty in colours. 

Plots of Wieros against the K magnitude for EROS LMC 
and SMC are presented in Fig. [9] Separation of C- and O-stars 
becomes particularly clear for the LMC. The location of three 
groups of stars with respect to the TRGB is visible in these 
plots. Simila r diagrams for the We senheit index were alr eady 
presented bv lNoda et"ai1 d2002l) and lSoszvnski et all (J2009), but 
our filter bands and fields are rather different. In particular the 
accuracy of the estimation of the Wesenheit index is sensitive to 
the quality of the two-band coverage. For EROS the advantage 
stems from the simultaneous observation in two bands. 

As the Wesenheit index is well correlated with K magnitude, 
it is in principle possible to implement the method of Sect. 14.31 
by relying purely on EROS visual observations. In particular the 
Wesenheit index would be helpful to separate the RGB and AGB 
stars; however, such a classification would be less precise than 
than the one based on K luminosity. The gap on TRGB, clearly 
visible in K distribution, is hard to find in Wi. We can also sep- 
arate RGB from AGB stars using the I-band. The TRGB appear 
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Fig. 11. Slope ay against V - R amplitude ampvR (top) and 
against 2MASS K magnitude, for MACHO LPV V and R pho- 
tometry of the LMC. The C-rich stars are marked with circles. 



at / = 14.54 in the LMC and / 
2000l). 



4.6. OGLE. 



14.95 in the SMC (ICioni et alJ 



The bands employed by the EROS survey were non-standard. 
In this respect it is desirable to verify our results with OGLE 
observations obtained in the standard / and V bands. For the 
OGLE light curves, we performed the calculations similar to 
those for EROS. In this way we obtained values of the corre- 
sponding slope ay. The results for OGLE are pres e nted i n Fig. 
[ITfl We only plot LPV selected bv lSoszvnski et all d2005l) . who 
reject low-amplitude variables. Thus the number of stars with 
low colour amplitude and low K luminosity is lower than for 
EROS LMC. 

Nevertheless, features revealed in Fig. |4]are present in Fig. 
[TOl It is remarkable that, despite the small number of V obser- 
vations and the interpolation of / magnitudes, separation of stars 
int o respective classes is c lear. 

Soszvnski et alJ d2005l) suggest that this sample should only 
contain the LPV AGB stars. We found that below the TRGB 
their stars follow the relation corresponding to the RGB stars. 
This may indicate that these stars belong to the AGB or at least 
that they follow the same relation as the RGB. These authors 
introduced the Period- Wesenheit index plots for a new method 
to differentiate between the O- and C-rich stars. Our method in- 



troduced in Sect. 3] is complementary as we use no information 
on periods. In this way our method is less demanding on time 
coverage of the light curves. Nevertheless, the two methods of 
selection of O- and C-rich stars are fairly consistent. 

4.7. MACHO. 

The same procedure as for EROS was repeated for the MACHO 
data. MACHO uses non-standard filters fairly close to the stan- 
dard V and R rather than V and /. In this way MACHO bands 
differ markedly from those of EROS and OGLE. This difference 
strongly affects our results for MACHO. The plots of the slope 
ay against either the colour amplitude for MACHO bands or the 
2MASS K mag reveal no clear separation of the O- and C-rich 
sta rs in Fig.fTT] This r esult seems consistent with the failure of 
the IWood et all d2004l) classification of the O- and C-rich stars 
based on similar filters. Thus our method only works for suit- 
able photometric bands and is useless for MACHO photometry. 

4.8. Maps of CIO ratio for Magellanic Clouds. 

In Table[T]we compared the three new methods described in Sect. 
|4T]|43] and [44] with the well known method K(J - K) that uses 
J — K > 1 .4 as a criterion for C-rich stars. Values in Table [T]were 
calculated only for variable stars. For all methods and for both 
Magellanic Clouds we found more C-rich stars than with the old 
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Fig. 10. Separation of C-rich, O-rich AGB, and red variables below the TRGB for LMC from OGLE LPV I and V photometry 
(ISoszvnski et al.ll2005l) . The colour slope parameter ay vs. the amplitude of colour change (left panel). The colour slope parameter 
ay vs. K 2MASS magnitude (middle panel). Colour slope parameter vs. J - K colour (right panel). 



Table 1. The C/O ratio for the MC's from variable stars 



Method 


K( J - K) 


a v (amp vl ) 


a v (K) 


a v (J - K) 


LMC 


0.63 


0.69 


0.69 


0.72 


SMC 


0.46 


0.89 


0.98 


0.81 



method. The K( J-K) method yields the conclusion that the C/O 
ratio is much lower for the SMC than for the LMC. The a v ( J-K) 
method yields the opposite result: the C/O ratio is greater in the 
SMC than in the LMC. Right panels of Fig. [4] shows that many 
C-rich stars lie between 1.15 and 1.4 on the J—K axis, especially 
for the LMC and a y(J - K) can separate stars without problem in 
this range of J - K. Different distributions of C-rich stars in the 
Magellanic Clouds probably depend on metallicity. Stars in the 
SMC are bluer than in the LMC for a similar stage of evolution. 

Table[T]also contains results for av(ampvi) and ay{K) meth- 
ods. To compare ay{ampvi) with other methods, we divided O- 
rich stars to RGB and AGB stars using TRGB. Both methods are 
affected by overlapping O- and C-rich regions but still we found 
more C-rich stars than with the the K(J - K) method. 

From our samples of O- and C-rich AGB stars from the 
ay(J - K) method, we calculated the C/O ratio in all cells be- 
longing to a grid of 100' x 100' over the face of the LMC and 
grid of 45' x 45' over the face of the SMC. Figure fT2lshows the 
ratios after boxcar average smoothing. A large hole in the central 
part of the SMC results from the lack o f data in that sector. We 
retreive the same trends as described by Cioni &Habind (120031) 
in their Figs. 3 & 4: in the LMC the higher values of the C/O 
ratio are located in the outer regions, while the structure in the 
SMC is clumpier. 
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5. Clues to the underlying scenarios 

Detailed modelling of the observed effect is beyond the scope of 
the present paper. As we pointed out above, modelling the late 
evolutionary stages is difficult and constitutes a whole branch 
of the theory. Worse, reconstruction of the synthetic spectra 
of C-rich AGB stars proved difficult an d strongly relies o n 
progress in the low-temperature opacities (Marigo et al. 2008). 
Simulation of the pulsations of the AGB stars is still more com- 
plex. The stellar at mosphere during pulse s does not move uni- 
formly up or down (llreland et al.ll2004blal) . Thus the non-linear 
pulsation model of a red giant is required to fit observatio ns 
(iLebzelter & Woodll2005l lOlivier & Woodll2005t rWoodll2006l) . 

However, it is possible that the effect demonstrated in Sect. 
|4] does not depend strongly on the details of the underlying 
physics. To demonstrate that we resorted to a crude and pos- 



Fig. 12. Maps of the C/O ratio in the LMC (upper panel) and 
in the SMC (lower panel) obtained from EROS data using the 
ay(J - K) method described in Sect. 14.41 



sibly non-physical model assuming that for all AGB stars dur- 
ing pulses (ilogL/iilog T e ff = const and that the amplitude of 
variation in L and T e fj remains the same for all stars. Thus any 
observed differences between stars would follow from a differ- 
ent response by their envelope to the same modulation. Next, we 
assume that the response of the upper envelope is quasi-static 
and may be determined by interpolation between different stat- 
ics, i.e. non-pulsating, evolutionary models from different tracks 
but with the suitable L and T e ff- This could be far from reality. 
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Fig. 13. Same as in Fig.|4]for the model slope parameter a v calculated forreal LMC stars using Marigo etall (120031) model sequences 
(see text for details). Our crude model yields a clear separation of O-rich RGB, O-rich AGB, and C-rich AGB stars. 



We employed the tracks calculated by iGirardi & Marigol 
(2007) and available for downloacQ to interpolate properties 
of individual stars at their extreme luminosities, we take lu- 
minosity L, effective temperature T e ff, absolute bolometric 
magnitude MboU absolute magnitudes in Johnson-Cousin-Glass 
UBVRIJHGK pass-bands, and surface carbon-to-oxygen num- 
ber ratio C/O listed for the tracks. We interpolated these values 
to find how luminosities in filters vary with L and T e /f during 
pulsation. In each track we identified the location of the RGB 
stars during TRGB and AGB phases of evolution. Next we per- 
turbed their L and T e ff to simulate the pulsation extreme phases. 
Using the corresponding filter magnitudes, we calculated slopes 
and amplitudes and plotted them in a similar way as for Fig. |4] 
By trial and error we found that a reasonable match with Fig. 
His obtained for AlogT e// = 0.042, AlogL = 0.026, which 
leads to dlogL/dlogT e ff = 0.614, but this works well for a 
wide range of both parameters. In this way we obtain Fig. [13] 
for the simulated data, where [13]0- and C-rich stars separate in 
the slope-colour amplitude, sloped luminosity and slope-colour 
J-K diagrams. 

As our model is entirely artificial, we refrain from strong 
conclusions. However, this result should encourage theoreticians 
to reproduce our e ffect using fully realistic calculations (e.g. 
Marigo et al. 2008), and to use it to restrict the pulsation sce- 
nario. 



6. Special properties of red variables. 

So far, our analysis has concerned the red variable stars with a 
modest amplitude (ampv < 1 mag.), since they posed biggest 
classification problems. Also we payed no special attention to 
the multiperiodic variables in general and those with LSP in par- 
ticular. Although our method for these stars does not explicitly 
account for their peculiarites, we briefly present results of its 
application hoping to provide additional clues missed in more 
specialised analyses. 
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Fig. 14. Colour-slope parameter in function of amplitude of 
colour change for EROS LMC and SMC data. Crosses mark 
large-amplitude variables. 



6.1. Large-amplitude pulsations 

In Fig. [14] we plot the slope- against colour amplitude for the 
whole range of amplitudes. This figure, which based on the raw 
data, corresponds exactly to the left column in Fig. |4] except for 
the extended amplitude scale. To illustrate here and in the next 
figure the relation of colour and luminosity amplitudes, ampvi 
and V, stars to the right of the line of inclination 0.6 passing 
through the origin are marked with different symbols. It appears 
that the pattern of the two populations of stars with small and 

4 http://pleiadi.aopd.inaf.it 



large slopes discussed in Sect. [4]does extend to the large ampli- 
tudes, too. On a related plot of Wesenheit index vs. K luminosity, 
the large-amplitude C-rich stars appear redder, while the O-rich 
ones are bluer then the corresponding stars of small amplitude. 
We found no other peculiarities for the large-amplitude stars. 

6.2. Long secondary pulsations 

Up to now we have analysed the raw data. However, many RGB 
variables reveal combinations of two periodic variations: a short 
one, presumably a radial pulsation, and a LSR The correspond- 
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Fig. 15. Correlation between amplitudes of LSP and of shorter Fig. 16. Correlation between the colour slope parameter of LSPs 
period variation for EROS LMC and SMC data. Crosses mark and the shorter period for EROS LMC and SMC data, 
large-amplitude variables. 



ing periods form sequences B/B' and D in the period- luminosity 
diagram (P-L). The natu re of the LSP in RGB stars is still a mys- 
tery for interpretation (ISoszvriskil (120071) : iNicholls etafl {2009) 
and references therein). The D sequence periods are roughly ten 
times longer than the B ones and correspond to no known ra- 
dial pulsation mode. The B sequence may correspond to a low- 
number overtone radial mode. LSP is observed among at least 
25% of low mass RGB stars and typical red amplitudes amp^ do 
not exceed 0.8 magnitudes. 

In this section we apply methods of Sect. [4] to the time- 
filtered light curves of all EROS LMC and SMC red variable 
stars, whether known to display an LSP or not. In the process 
we use the fact that our procedure does not require a known pe- 
riod. By employing the filters described in Sect. |3]we extracted 
separate light curves for fast and slow variations. The cut-off fre- 
quency of our filters depends on the mean magnitude to ensure a 
proper split of B and D sequences in systems known to display 
LSPs. Thus we obtain two different sets of amplitude and slope 



parameters ampv- s hort and ampv-i on g, av-shon and ay-long from 
respectively high- and low-pass filtered light curves. 

In Fig. [15] we plot long-time amplitudes of variation of mag- 
nitude against short time ones, ampv-i on g v.s. ampv-short for all 
EROS LMC and SMC data. Stars marked in Fig. OH] as pos- 
sessing relatively large colour magnitude ampvi also exhibit 
absolutely large magnitudes ampy on at least one time scale. 
Although the amplitudes in Fig. [15] were estimated indepen- 
dently from two different light curves, they seem correlated. A 
consistent yet non-unique explanation would involve a similar 
mechanism for the short and long light variations. 

A stronger argument supporting the similarity of mecha- 
nisms of short and long time variations stems from the tight cor- 
relation of the short and long time scale slopes a s h or t and ai ong re- 
vealed in Fig. [16] As discussed in Sect.[4]the two separate blobs 
in Fig.[T6lcorrespond to C- and O- rich stars. The presence of just 
two such blobs rather than three or four at square vortices pro- 
vides striking demonstration that short and long time scale vari- 
ations both depend in the same way on surface chemical prop- 



Mariusz Wisniewski et al.: Populations of variable red giants in the Magellanic Clouds 



13 



erties. Combined evidence from Figs. 1X31 and [TBI indicates that 
both long and short time scale variations are due to the reaction 
of the stellar photosphere to similar kind of perturbations . We 
caution to not overinterpret this result. Firstly, our data contain 
both stars with and without known LSPs, and most light curves 
exhibit both periodic and irregular variations. However, it must 
be remembered that LSP amplitudes are not that small, and they 
are present in a sizable fraction of stars. If the LSP were pro- 
duced away from the photosphere, one could expect to see some 
substructure within both C and O blobs in Fig. [16] so apparent 
lack of any internal structure in these blobs is consistent with 
LSPs being caused by the same kind of photosphere response as 
for the irregular variations and for the short period pulsations. 
In this context we stress that the EROS two-band observations 
are strictly simultaneous; hence our conclusions on photosphere 
evolution involve no interpolation or implicit light curve models. 



7. Conclusions. 

A variable star during each cycle completes a closed loop on the 
colour-magnitude diagram. Non-periodic variables travel more 
involved routes. For each red variable star from the EROS sur- 
vey, we derived a width and average slope of these figures, re- 
spectively ampvi and ay. In doing so, we benefitted from re- 
liable colours obtained from EROS simultaneous two-band ob- 
servations. In this way our analysis is less affected by variability 
and/or time sampling patterns. In particular, the gaps and alias- 
ing pose no big problem for us as long as the data cover most of 
the pulsation period. This is a much less demanding requirement 
than a span of several cycles needed for reliable period analysis. 
The price paid in the process is a loss of any extra information; 
nevertheless, our analysis yields independent results that supple- 
ment more traditional period studies. 

Our key method derived in Sect. 14. 1 I relies on the parameters 
ay and ampvi derived from pure visual observations. For its test- 
ing and verification we cross-referenced stellar positions from 
the EROS, OGLE and MACHO optical surveys with the infrared 
2MASS catalogue. K magnitudes from the last catalogue were 
treated as indicators of stellar mean luminosity. Such a proce- 
dure is reliable for a known distance because it has little K band 
extinction and small bolometric correction for red variables. To 
completely free our optical luminosities from extinction, we em- 
ployed the Wesenheit index. However, our key parameters result 
ay and ampyj are independent of distance and magnitude val- 
ues. The range and correlation of two band magnitude changes 
are the only parameters that are important for us. 

We demonstrated that in the ay and ampyj diagram, red vari- 
able stars form two well defined bands along ay = 0.4 and 
0.65 lines. The bands occupied the same location for both the 
LMC and SMC stars. We argued that such a diagram consti- 
tutes a new method for differentiating of C- and O-rich RGB 
stars from purely visual observations. So far, effective photomet- 
ric methods have relied on infrared observations. In fact, the di- 
chotomy of C/O star distribution becomes even more prominent 
in the slope ay- K luminosity and slope- (J-K) colour diagrams 
for all stars from our sample. The clustering is particularly pro- 
nounced, both for LMC and SMC, in the slope-luminosity dia- 
gram. There the AGB and RGB stars separate well, apart from 
the usual C/O split of ay. In this particular diagram stellar evo- 
lutionary tracks RGB -TRGB- AGB (O)- AGB (C) are particularly 
well defined. While the detailed location of clusters in our di- 
agnostic diagrams does depend on passbands, it is clearly man- 
ifested both for both EROS and OGLE two-band photometries. 



However, the bands employed by MACHO are not suitable to 
our classification method. 

The final proof of the physical meaning of the clustering 
observed in our diagrams stems from the cross-check with cat- 
alogues of C-rich stars. The overwhelm i ng ma j ority of C-rich 
stars fr om catalogues by iKontizas et al.l (l200lh. iGroenewegenl 
(l2004t) . lRebeirot et all dl993l) . lMorgan & Ha tzidimitrioul (ll995l) . 
for LMC and SMC, belongs to the low- ay slope cluster. Thus 
it appears that indeed the C- and O-rich stars belong to separate 
clusters in our ay vs. ampyj diagnostic diagram. In this way we 
discovered and verified a new method of distinguishing between 
C- and O-rich stars. Our method works for sparse coverage in 
two-colour light curves in the bands resembling V and /. Since 
we rely on colours, the measurements in both filters should be 
(nearly) simultaneous. Our method was tested successfully for 
EROS-2 LMC and SMC observations and also for OGLE data. 
It does not seem to suffer from any metallicity effect as clus- 
ter locations for the LMC and SMC were identical. Using sim- 
ple simulations based on realistic evolutionary models of RGB, 
we demonstrated plausible cause of the effect revealed by our 
method; however, for as broad filters as those used by MACHO 
our method fails. One future application of our diagrams is to 
test the pulsation evolutionary models on a large population of 
RGB and AGB stars. 

Using our method of differentiating the O- and C-rich stars 
we were able to derive the population C/O star ratio, an indica- 
tor of the mean metallicity of a stellar population, as wel l as a 
tracer of the history of stellar formation (Cioni et al . 2006). Our 
selection can recognise more C-rich stars than using J - K > 1 .4 
criterion, and it yields higher values of C/O ratio. Thus we were 
able to produce C/O maps for the LMC and SMC. 

We payed some attention to colour correlation properties of 
the LSP compared to the short primary pulsation (SPP). Our re- 
sults provide some evidence favouring origin of both LSP and 
SPP in the stellar photosphere. This follows from the correla- 
tion of colour-magnitude slopes and separately from correlation 
of amplitudes for the LSP and SPP. In other words, both SPP 
and LSP depend on the chemical (C/O) properties of the pho- 
tosphere. The colour variations would be consistent with mod- 
ulation of the effective temperature. This seems to exclude any 
aspect effects, e.g. eclipses, as a cause of LSP; however, using 
our purely photometric observations it would be premature to 
conclude that both SPP and LSP clocks are due to pulsation. 
However, both clocks should somehow affect the photosphere in 
similar ways. 
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